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Abstract 

We report on the results of scanning tunneling spectroscopy measurements on single 
crystals of MgB2. Tunneling was performed both parallel and perpendicular to the 
crystalline c-axis. In the first case, a single superconducting gap (At^ = 2.2 meV) 
associated with the vr-band is observed. Tunneling parallel to the a6-plane reveals an 
additional, larger gap (Ao- ~ 7 meV) originating in the highly two-dimensional a- 
band. Vortex imaging in the vr-band was performed with the field and tunnel current 
parallel to the c-axis. The vortices have a large core size compared to estimates 
based on -ffc2) and show an absence of localized states in the core. Furthermore, 
superconductivity between the vortices is rapidly suppressed by an applied field. A 
comparison to specific heat measurements is performed. 
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1 Introduction 



Since the first report of superconductivity with a remarkably high Tc = 39 K 
in magnesium diboride (MgB2) only two years ago [1], enormous progress has 
been made in the studies of this material. From a fundamental point of view, 
a central issue has been to establish that MgB2 is a two-gap superconductor 
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- a concept which was introduced already in the fifties [2,3], and which has 
now found renewed relevance. Two-gap superconductivity in MgB2 was first 
predicted theoretically [4] and is now commonly accepted, confirmed by a large 
number of experiments [5,6,7,8,9,10,11,12]. 

The two gaps in MgB2 are associated with different parts of the Fermi surface, 
which is composed of 4 seperate sheets [4,13]. Two of these are coaxial cylin- 
drical sheets parallel to c*, derived from cr- antibonding states of the boron p^y 
orbitals (cr-band). The other two sheets are derived from 7r-bonding and an- 
tibonding states of the boron orbitals (7r-band) and are three-dimensional. 
Calculations of the superconducting gap size for the two different bands yields 
Ao- 7 meV, and = 1 — 3 meV [14,15]. These values are respectively 
above and below what is expected for a BCS s-wave superconductor, A ~ 
1.76 ksTc — 5.9 meV. Furthermore, and in contrast to many materials or al- 
loys studied earher, the two bands in MgB2 have roughly equal weight, leading 
to new and interesting phenomena as we will show in the following. In this re- 
spect it has become clear, that local spectroscopic investigations of the mixed 
state is an ideal way to study the detailed nature of the superconducting state 
in MgB2 at a microscopic level. 

In this paper we summarize the present status of our scanning tunneling spec- 
troscopy measurements, made possible by the recent availability of high qual- 
ity MgB2 single crystals. The results include spectra obtained by tunneling 
parallel as well as perpendicular to the crystalline c-axis, showing how the 
tunnel direction changes the coupling to the two bands. In addition, vortex 
imaging was performed with the tunnel current and magnetic field parallel to 
the c-axis. The vortices are found to have a number of remarkable properties: 
An absence of localized states, a very large vortex core size compared to the 
estimate based on ifc2, and a strong core overlap. Some of these results were 
reported earlier [16]. 



2 Experimental 

The scanning tunneling spectroscopy (STS) experiments were performed using 
a home built scanning tunneling microscope (STM) installed in a ^He, UHV 
cryostat containing a 14 T magnet [17]. The tunneling measurements were 
done using electrochemically etched iridium tips, with the bias voltage applied 
to the sample, i.e. a positive bias corresponds to probing the empty states and 
negative bias to the occupied states respectively above and below the Fermi 
level. A tunnel resistance of Rt = 0.2 — 0.4 Gf2 was used, and the differential 
conductance measured using a standard ac lock-in technique. 

Single crystals of MgB2 were grown using a high pressure method described 
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elsewhere [18], yielding platelike samples with the surface normal parallel to 
the crystalline c-axis. The surface area of the crystals are roughly 0.25 x 0.25 
mm^, and the thickness varies between a few up to ~ 100 microns. The critical 
temperature is typically Tc = 37 — 38.6 K, with a sharp transition, AT^ = 
0.4 — 0.6 K, measured by SQUID magnetometry [18]. Measurements were done 
both with the tunnel current. It, parallel and perpendicular to the c-axis. For 
It II c, the spectroscopy was done on the surface of an as grown single crystal. 
In the case of /t -L c, a relatively thick crystal was cracked to expose a clean 
surface and immediately mounted and inserted into UHV. 

The upper critical field of MgB2 is highly anisotropic [19], with zero temper- 
ature extrapolations of respectively 11^2 = 3.1 T and i/^^ = 23 T for single 
crystals [20], even though the latter value is subject to some uncertainty. Using 
the Ginzburg-Landau (GL) expression for ifc2 = 0o/(27r^^), where 0o = h/2e 
is the flux quantum, yields an in-pane coherence length, — 10 i^™- to 
the large anisotropy, .^ql — 'Cgl ^ ^02/ ^c2 — 1 — 2 nm. An estimate of the 
mean free path, based on the measured residual resistivity [21] and specific 
heat [11], and the calculated Fermi velocity [22], gives / = 50 — 100 nm, and a 
flrst estimate is thus that the samples are in the moderately clean limit [23]. 
However, the two-band nature of MgB2 comphcates the picture since IIc2 is 
dominated by the two-dimensional a-band which gives rise to the large ifc2- 
anisotropy, while the transport properties depends on both bands. A detailed 
analysis of the thermal conductivity by Solugubenko et al. concluded that the 
mean free path is roughly equal for the two bands with I ^ 80 nm [21]. 



3 Results 

In the following we present the results of our measurements on MgB2 single 
crystals. 

3.1 Dependence on the tunnel current direction 

The observation of a single or both superconducting gaps depends on the ori- 
entation of the sample, since tunneling along different directions changes the 
coupling to the two bands [22]. Roughly speaking, only bands with a compo- 
nent of the Fermi velocity parallel to the tunneling direction are observed. In 
particular, this means that coupling to the 2D cr-band is highly suppressed for 
tunneling parallel to the c-axis, whereas the 3D 7r-band will be probed for all 
tunneling directions. 

In Fig. 1(a) we show a superconducting spectrum with A || c. This is an average 
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Fig. 1. Zero- field superconducting spectra of MgB2 obtained by tunneling respec- 
tively parallel (a) and perpendicular (b) to the c-axis. In both cases T = 0.3 K and 
Rt = 0.2 GQ (U = 0.1 V; / = 0.5 nA). 



of 40 spectra obtained along a 100 nm path, which shows perfect homogeneity. 
One observes a single gap with coherence peaks at ±2.9 meV associated with 
the TT-band, and additional weak shoulders at ±6 meV as indicated by the 
arrows. The shoulders are an admixture of the gap in the cr-band. 



In Fig. 1(b) we show an average spectrum obtained for It _L c. For this di- 
rection of tunneling, both superconducting gaps are present. The position of 
the low-energy coherence peak is unchanged as shown by the dashed lines. 
The additional peak arising from the cr-band is located at ±7.3 meV. This 
is the first report of tunneling measurements on single crystals of MgB2 for 
both /t II c and It ± c, which allows a direct correlation between the tunneling 
direction and the observed gap(s). Our results agree with those of lavarone et 
ai, who investigated a number of single grains with different, but unknown 
absolute orientations [12]. 

For both tunneling directions the low zero bias conductance indicates a high 
quality tunnel junction and a low noise level. True vacuum tunneling condi- 
tions were confirmed by varying the tunnel resistance, Rt, and verifying that 
the spectra normalized to the conductance outside the superconducting gap 
collapse onto a single curve. 
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3.2 Temperature dependece for It\\c 



Tunneling parallel to the c-axis, we measured the superconducting spectrum 
for a number of temperatures between 320 mK and 38.8 K. As discussed 
above this probes mainly the 7r-band. The results are presented in Fig. 2(a), 
and shows how the gap is gradually suppressed and finally closes at Tc. We 
have fitted the measured spectra to the conductance calculated using the BCS 
expression for the density of states, including a finite quasiparticle lifetime, F 
[24]: 

e-ir 



N{e,T)=No 
The conductance is given by 



Re 



(1) 



.,.)^^./.e.(e.r)(-ML±£lO). 

which we furthermore have convoluted by a Gaussian of width u to take into 
account experimental broadening. The fits are shown in Fig. 2(a) and are found 
to be in good agreement with the measurements, except at low temperatures 
where the height of the coherence peaks is overestimated. The extracted values 
of the gap are shown in Fig. 2(b), yielding A7r(0) = 2.2 — 2.3 meV. We find 
that A^(T) is in excellent agreement with the temperature dependence of 
the BCS A(T) [25] when the zero temperature gap is scaled to 2.3 meV. The 
lifetime parameter F = 0.1 meV at low-T and increases with temperature. The 
experimental broadening u; — 0.5 meV at base-T and increases only slightly to 
0.8 meV at 10 K. For comparison the ac excitation used for the measurements 
was 0.4 meV RMS. Above 10 K the effects of F and u; becomes identical, and 
the latter is kept fixed. 



3.3 Vortex imaging in the n-band 



We now turn to measurements in an applied magnetic field. These measure- 
ments were done with It \\ H \\ c, and hence selectively probes the 7r-band. 
In a type-II superconductor such as MgB2, the magnetic field penetrates into 
the sample in the form of vortices each carrying one fiux quantum, and gen- 
erally arranged in a periodic array: the vortex lattice. In the core of each 
vortex, superconductivity is suppressed within a radius roughly given by the 
coherence length, ^. The vortex spacing is determined by the applied field 
and the fiux quantization, and in the case of a hexagonal vortex lattice it is 
d^{2/^3cf>o/Hy/\ 

Necessary prerequisites for obtaining vortex lattice images is a high sample 
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Fig. 2. (a) Spectra obtained a number of temperatures between 320 mK and 38.8 
K (thick lines), tunneling parallel to the c-axis. Each spectrum is normalized to 
the conductance at 12 meV, and offset by 0.5 with respect to the previous one 
for clarity. The thin lines are our fits to the spectra as described in the text, (b) 
Temperature dependence of the superconducting gap in the vr-band. The line is the 
BCS expression for A(0) = 2.3 meV and Tc = 39 K [25]. 

homogeneity and a relatively flat sample surface. We verified that the first 
requirement is fulfilled by measuring perfectly reproducible spectra in zero 
field at a large number of positions along a several thousand angstr0m line. A 
topographic STM image of the sample surface is shown in Fig. 3. The surface 
is seen to be very smooth, with an indication of a granular structure with a 
characteristic size of 5 - 10 nm. No effect of this "granularity" is observed in 
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100 nm 



Fig. 3. STM topographic image of the as-grown surface of a MgB2 crystal. The 
RMS roughness is 3 A. 




250 nm 250 nm 200 nm 



Fig. 4. (color) Vortices in MgB2 at respectively 0.05 T(a), 0.2 T(b) and 0.5 T(c) at 
2 K. The spatial variation of the conductance is shown by the color scale (different 
for each image). 

the spectra. 

The magnetic fields were apphed at 2 K, and the system allowed to stabilize for 
at least a few hours. After this time no vortex motion was observed, indicating 
a fast relaxation and hence low vortex pinning in the crystal. In Fig. 4 we 
show STS images of vortices induced by three different applied fields. The 
images were obtained by measuring the differential conductance at zero bias. 
Low values of the conductance correpond to superconducting areas, and high 
values to the vortex cores. Figs. 4(b) and (c) show a well ordered hexagonal 
vortex lattice, with a lattice constant that corresponds to the applied field 
within 10%. 

The low field in Fig. 4(a) is equivalent to a separation, d = 220 nm ^ ,^gl- The 
vortices can therefore be considered as isolated from each other. Such isolated 
vortices are expected to contain localized quasiparticle states (Andreev bound 
states), which should show up as a zero bias conductance (ZBC) peak at the 
vortex centre [26,27], provided that the sample is sufficiently clean to prevent 
these from being smeared out by scattering. We have measured the evolution 
of the spectra at a large number of positions along a trace across the vortex 



7 



^to ^5 6 5 10 
Bias voltage (mV) 




n.n l 

-150 -100 -50 50 100 150 

Position (nm) 

Fig. 5. (color) (a) 250 nm trace through the centre of the single vortex in Fig. 4(a), 
with spectra recorded each 2 nm. A spectrum at the vortex centre together with 
one far from the vortex core have been highlighted in red for clarity, (b) Normalized 
zero bias conductance versus distance from the centre. The line is a fit to eq. (3) in 
the text. 

core. This is shown in Fig. 5(a). Contrary to expectations, we find that the 
normalized ZBC increases to one with no indication of any localized states. 
Instead, the spectra at the centre of the vortex are absolutely flat, with no 
excess spectral weight at or close to zero bias. This absence of localized states 
is striking, considering that / = 5 — 10 x ,^gl- However, as we will show below, 
the coherence length in the vr-band is approximately 50 nm. This is much larger 
than the GL estimate based on Hc2, and equal to only 1-2 times the mean 
free path. Nonetheless, for comparison systematic studies of Nbi_a;Ta2;Se2 with 
X = — 0.2, showed that even for ^/Z 1 some excess weight close to zero 
bias was observed [27]. 

The spatial extent of the vortex is shown in Fig. 5(b), where we have plotted 
the normalized ZBC, cr'(x, 0) from the vortex trace. It is immediately clear 
that the spatial extension of the vortex core is much larger than the 10 nm 
estimated from Hc2 ■ We find that the vortex profile can be fitted by one minus 
the GL expression for the superconducting order parameter: 



a'{x,0) = o-Q + (1 - o-q) X (1 - tanhx/0, 



(3) 
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Fig. 6. Normalized spectra at 2K measured in zero field, and between the vortices 
for fields between 0.05 T and 0.5 T. Each subsequent spectrum is offset by 0.5 with 
respect to the previous one, with the bars at zero bias indicating the respective zero 
conductance. 

where ctq = 0.068 is the normalized ZBC measured in zero field. The fit yields a 
coherence length, ^ = = 49.6 ±0.9 nm. Using the GL expression to calculate 
the upper critical field with this value of the coherence length, yields H^2 — 
0.13 T. This value is in strong contrast to the fact that we have performed 
vortex lattice imaging up to 0.5 T, and shows that a separate Hc2 for the 
TT-band does not exist. 

The large has however other consequences. In Fig. 6 we have plotted the 
normalized ZBC between the vortices. From this it is clear that even modest 
fields rapidly suppress superconductivity in the region between the vortices, 
seen both by an increase of the ZBC and by a suppression of the coherence 
peaks outside the vortex cores. This is consistent with earlier point contact 
spectroscopy measurements [6], with the addition that we resolve the local 
behaviour on a microscopic scale. Such a behaviour is only expected for fields 
close to Hc2, and corresponds to a significant core overlap. 



4 Discussion 

The vanishing of at the bulk Tc and the sustained superconductivity de- 
spite a large vortex core overlap for H > 0.2 T leads to the conclusion that 
superconductivity in the vr-band is induced by the cx-band, either by interband 
scattering, or Cooper pair tunneling [3,28]. This means that by itself the vr- 
band would either be non-superconducting or have a very low and upper 
critical field. Consequently the observed behaviour reflects the state in the 
(T-band by an interband proximity effect, along the lines of recent theoretical 
work [28] . This is also consistent with estimates of the coherence lengths, using 
the BCS expression = hvp/^TiA^O)) and considering each band separately. 
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Taking the calculated average Fermi velocity in the a6-plane for the 7r-band. 
Vp — 5.35 X 10^ m/s [22], and the measured gap value A^(0) = 2.2 meV we get 
^0 — 51 nm, in excellent agreement with obtained from the vortex profile. A 
similar analysis for the cr-band, using Vp = 4.4 x 10^ m/s [22] and A^{0) ~ 7.1 
meV [12] yields — 13 nm. This agrees with the coherence length obtained 
from Hc2, and reinforces the conchision that it is mainly the cr-band which is 
responsible for superconductivity in MgB2, and thus determines the macro- 
scopic parameters T^. and Hc2- One must however bear in mind that this is a 
very simple-minded analysis. A detailed theoretical treatment of two weakly 
coupled superconducting bands in the presence of vortices, and the possibility 
of having have different coherence lengths, have to our knowledge not yet been 
presented. 

4-1 Comparison to specific heat measurements 

The vortex core overlap explains the strongly nonlinear field dependence of the 
electronic specific heat (T-linear term), 7 [5,11]. In type-II superconductors 
core overlap is usually negligible, and in the simplest approximation each 
vortex creates the same number of quasiparticles at the Fermi surface, hence 
contributing by the same amount to the specific heat. In that case 7 = 7n x 
H/Hc2, where 7n is the electronic specific heat in the normal state. However, in 
the case of MgB2, with strong core overlap in the 7r-band, the isolated vortex 
assumption is violated. Constructing a very simple model for the core overlap 

by 

a'(r, 0) = a(, + (1 - a'^) x (^1 - n,tanh , (4) 

where Vi are the vortex positions for a hexagonal lattice with a density corre- 
sponding to the applied field, we can calculate the ZBC at any position in the 
vortex lattice unit cell. In Fig. 7 we compare the calculated conductance at the 
midpoint between three vortices with the measured "bulk" ZBC. This shows 
a very good agreement, especially taking into account that there are no free 
parameters in the calculation: is determined from the vortex profile, and cTq 
is given by the finite ZBC in the zero field measurement. The 7r-band contri- 
bution to the specific heat can now be calculated by averaging the normalized 
ZBC in one vortex lattice unit cell, 7,^ = 7^ 0)) [28]. For this calculation 

we have set o"q equal to zero. On the other hand, the cr-band can be described 
by the usual linear field dependence 7cr = 7n x H/Hc2- Adding the terms gives 
7 = 77r -I- 7(T, where 7n/7n is the relative weigth of the two bands. The calcu- 
lated field dependence of 7 is shown in Fig. 7, in perfect agreement with the 
measured specific heat for polycrystalline MgB2[ll], using 7^/7^ = 0.55/0.45. 
The distribution of weight between the two bands is in agreement with other 
theoretical and experimental estimates [4,9,11,15,21]. Furthermore, the suc- 
cessful correlation with specific heat measurements makes us confident that 
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Fig. 7. Calculated "bulk" ZBC (left axis) and electronic specific heat, 7, for 
7n/7n = 0.55/0.45 (right axis). The calculated values are compared to respectively 
the measured bulk ZBC (circles), and specific heat measurements on polycrystalline 
samples (squares) [11]. 

our results reflect bulk properties of MgB2. 



5 Summary 

In summary, we have presented STS data on the vr-band in MgB2, including 
the flrst vortex imaging in this material. We have demonstrated the absence 
of localized states in the vortex core, a very large vortex core size and a strong 
core overlap. The data presents a striking experimental demonstration of the 
fundamentally different microscopic properties of the two bands in MgB2. 
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